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According to classical grain growth laws, grain growth is driven by the minimization of surface 
energy and will continue until a single grain prevails. These laws do not take into account the 
lattice anisotropy and the details of the microscopic rearrangement of mass between grains. Here 
we consider coarsening of body-centered cubic polycrystalline materials in three dimensions using 
the phase field crystal model. We observe as function of the quenching depth, a cross over between 
a state where grain rotation halts and the growth stagnates and a state where grains coarsen rapidly 
by coalescence through rotation and alignment of the lattices of neighboring grains. We show that 
the grain rotation per volume change of a grain follows a power law with an exponent of —1.25. 
The scaling exponent is consistent with theoretical considerations based on the conservation of 
dislocations. 


INTRODUCTION 

In polycrystalline materials, the microstructure, given 
by the grain sizes, shapes and coordination, is an impor¬ 
tant ’’control parameter” for macroscopic material prop¬ 
erties like the yield stress, conductivity and brittleness. 
Polycrystalline materials typically form from the nucle- 
ation and growth of grains with different lattice orienta¬ 
tions in a quenched or annealed melt. If the annealing 
temperature is sufficiently high, the grain boundaries will 
be mobile and rearrange in a way that lead to an overall 
coarsening of the polycrystalline matrix with time. 

In classical models of grain growth, like the Neumann- 
Mullins model, grain boundary migration is driven by 
surface energy and the average grain area grows linearly 
with tim^, (A) ~ t. In three dimensions, however, a 
correpondingly simple universal growth law can not be 
establishecP^^ because of the increased geometrical com¬ 
plexity. For models whose dynamics are driven by the 
minimization of a free energy constrained with a con¬ 
tinuity equation, general arguments relying on energy 
dissipatioiP sets the upper bound for the coarsening dy¬ 
namics in D dimensions to {L^) ^ which coin¬ 

cides with the Neumann-Mullins law. Experimentally it 
has been observed that the grain size growth can be de¬ 
scribed by a power law with an exponent that depends 
on a variety of experimental factors such as the annealing 
temperatur^. It is further found that the grain growth 
is sublinear, i.e. the growth in area is described by a 
power law exponent less than unitjPl. Another possible 
mechanism, which results in grain coarsening, is the coa¬ 
lescence of neighbouring grains. The rotation of crystals 
might first lead to the lattice alignment of neighboring 
grains and secondly the elimination of the grain bound¬ 
aries between them. Highly convoluted grains can be 
formed in this way. Grain rotation has been observed 
in experimental and has recently also been observed in 
molecular-dynamics simulation^. 

In many materials, the coarsening dynamics further 
differs from normal grain growth, in the sense that the 
grain size distribution varies in time. These abnormal 
grain growth systems are characterized by a minority 


of ’’abnormal” grains growing faster than the mean size 
leading to an inhomogeneous size distributiorP. In con¬ 
trast to classical grain growth laws, grain growth are also 
known to stagnate in tim^^Sl, jn two dimensional sys¬ 
tems, it has been suggested that the stagnation might be 
a combination of high kinetic barriers, preventing mass 
migration across grain boundaries, and a locking of in¬ 
dividual grains preventing grain rotation and subsequent 
lattice alignmeniPl. 

In this article, we present a numerical study of the 
influence of grain rotation on the coarsening dynamics in 
three dimensional polycrystalline systems with a body- 
centered cubic (BCC) symmetry. 

MODEL AND ANALYSIS 

To study the coarsening dynamics, we use the phase 
field crystal model (PEG). The PEG model describes the 
evolution of a continuous order parameter field, which is 
spatially periodic with atomistic resolution. The model is 
based on the minimization of a phenomenological Swift- 
Hohenberg free energy functional given by 

-7^ = / dr Qv- (1 + V' + y + y ^ (1) 

with the order parameter if representing the crystal den¬ 
sity fielcP^. The parameter 02 and the mean density '0 
are related to the melting temperature according to the 
phase diagram in Fig. The points in Fig. [^show the 
parameters used in the simulations presented in this ar¬ 
ticle. While we change both the mean density and the 
parameter 02 in our simulations presented here, we shall 
for convenience only refer to the 02 value in the figures 
and in the text. The corresponding values for the mean 
density can then be found from the phase diagram. Al¬ 
ternatively the functional form of Eq. ( fll) c an be derived 
from classical density functional theorj®l. In three di¬ 
mensions a free energy on this form results in a rich phase 
diagram with multiple equilibrium phases. Here, we shall 
focus only on the BCC phase. The evolution equation of 
the density field 1 /) is assumed to obey an over damped 
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diffusion equation 


# _ V72 

dt 6^ ^ ^ 

Thereby establishing a link between the microscopic 
length scales and diffusive time scales. The PFC model 
has been used to study a wide range of phenomena in¬ 
cluding phase transitionJi^, plastic deformatiorfi^ and 
has been shown to successfully predict grain boundary 
energy as a function of mis-orientatioiP^l. We solve the 
dynamical equation using an exponential time integra¬ 
tion schem^i^ in a three dimensional box of size Lx LxL 
with L = 512dx or L = 1024dx. Time and space are dis¬ 
cretized with dx = ^ and dt = 1/2. We initialize our 
system from an undercooled melt by introducing small 
crystal seeds at random points and with random lat¬ 
tice orientations. After the initial crystallization phase, 
where the seeds grow to cover the whole melt, a poly¬ 
crystalline structure is formed which coarsen over time. 
During the coarsening stage, we track the volume and 
lattice orientation of all grains. 

For the segmentation of the grains, the peaks of the 
density field are located by a thresholding procedure from 
which the coordinates of the center of mass of each peak 
can be calculated. Grain boundary detection is per¬ 
formed using a Voronoi tessellation and a centro sym¬ 
metry parameteili^. The angle and axis of rotation for 
each grain is found by taking the mean over the orienta¬ 
tion of all the unit cells in the interior of the grain. The 
misorientation angle of two grains A and B is calculated 
by constructing the rotation matrix G for the grains and 
form the product AGab = GaG]^^. The misorientation 

angle 6 is then given hy 9 = arccos ^ 

To get a measure of the typical grain size in our 
systems, we calculate a coarsening parameter by fit¬ 
ting a Lorentzian squared to the structure factor of the 
systemP^, (|F['i/;]p). Where F[tp] is the Fourier transform 
and the outer brackets denote averaging over all orien¬ 
tations in fc-space. The width of the averaged 

spectrum provides a measure of the ordering scale of the 
system, which is inversely proportional to the mean grain 
size. 


RESULTS 

In general, the polycrystalline structure coarsen at a 
rate which depend on the quenching depth, the 02 pa¬ 
rameter in Eq. Q. In Fig.f^ we show the change in the 
ordering length, — (i, as a function of time 

for different quenching depths. Where is the initial 
ordering scale measured right after the system has fully 
crystallized. For deep quenching parameters, the system 
dynamics is described by a short period of grain growth 
followed by stagnation. As the quenching depth is de¬ 
creased, the time it takes to reach the stagnation stage is 
increased and eventually (for 02 ^ —0.15 and "(/j > —0.22 



FIG. 1. Mean size of grains on double logarithmic axes. The 
inset shows the same data on semi logarithmic axes. 


in the phase diagram. Fig. ^ the stagnation stage is re¬ 
placed by a stage of rapid grain growth. By tracking 
the orientation of individual grains, we observe that the 
stagnation is accompanied by a general decrease in the 
change of the lattice orientation. The change in orienta¬ 
tion eventually drops to zero when the system stagnates, 
see Fig. 
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FIG. 2. The PFG phase diagram with points indicating the 
parameters used in the simulations. The region in the middle 
has a BGG lattice as its equilibrium state while the equilib¬ 
rium states of the regions to the left and right are given by a 
uniform liquid state and a rod state, respectively. The gray 
areas indicates the coexistence regions. 

Independent of the system size, the average grain size, 
here represented by A^, can in the early state be ap¬ 
proximated by a power law in time, (A^) ~ with an 
exponent n that depends on the quenching parameters 
02 and tp. From simulations, we find values of the expo¬ 
nent in the interval n ~ 0.5 — 1. In order to estimate the 
variation of the coarsening exponent, we have performed 
repeated simulations for the same quenching depths and 
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FIG. 3. Average rotational change as function of time. The 
straight line corresponds to a power law with an exponent 
n = — 1. The stagnating regime shows a rapid decrease in 
grain rotation in contrast to the fast coarsening regime where 
grain rotation continues. 


we find a standard deviation in the power-law exponents 
of the order of 0.1. 

In general, the coarsening dynamics exhibits abnormal 
grain growth where a few grains grow significantly faster 
than the rest of the crystal matrix as can be seen in Fig.|^ 
The figure shows a late stage snapshot of a L = 1024:dx 
system. A consequence of abnormal growth is that the 
grain size distribution is not self-similar as can be seen 
from Fig. where the ratio of the standard deviation 
of the grain sizes and the mean grain volume is plotted. 
It is also clear from this figure that there is a strong 
dependence of the rotation on the parameter 02 . 

The transition to fast grain growth is initiated by the 
mobilization of small grains, which align in lattice ori¬ 
entation with the larger so-called abnormal grains and 
finally coalesce in to an even larger grain. This is sup¬ 
ported by the observation that the amount of rotation 
per change in volume V increases as the small grains get 
even smaller. From the scatter plots of |A0/AF| ver¬ 
sus F, we fit the general trend by a power law whose 
exponent is plotted in Fig. Best fits across multiple 
quenching depths and initial conditions suggests a scal¬ 
ing behaviour 


A6» 


y- 


( 3 ) 


where the scaling exponent is estimated to /3 = 1.25 ± 
0.06. The uncertainty is taken to be the standard devia¬ 
tion found by fitting a straight line using a weighted least 
squares approach. 


DISCUSSION 

In our simulations, we observe two distinctly different 
types of coarsening dynamics. For deep quenches the 



FIG. 4. Snapshot of a L = 1024da; simulation performed with 
02 = —0.15 and with periodic boundary conditions. Note the 
heterogeneous grain distribution which can be observed at the 
transition to abnormal grain growth. 


A 

g 

f-o 


Oa2 = -0.14 
□ a2 = -0.12 
0 a2 = -0.095 


□ 

^ O 


□ 

o 


□ 

.0 






o 

□ <p 


o 

o 


80 


□ o^ 




0.008 0.01 0.012 0.014 0.016 0.018 

(^) 


FIG. 5. The ratio of the standard deviation in grain size and 
the mean size as a function of mean size. The non-linear 
relation is a sign of abnormal grain growth. 


grain growth is typically accompanied by a decreasing 
rotation of the individual grains. When the grain rota¬ 
tion stops, the system has reached a stagnation state. 
Similar resu lts have been reported in two dimensional 
system^Ali8|_ shallow quenches, we find that the early 
stage coarsening is described by a power law growth of 
the form V oc Af ~ with n ~ 0.5 — 1, in reason¬ 
ably agreement with experimental studied, which finds 
n ^ 0.4 — 1.2. Further the grain growth is found to be 
abnormal resulting in a crystal matrix consisting of a 
few large convoluted grains. Abnormal grain growth has 
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FIG. 6. Powerlaw coefficients obtained from best fits to scat¬ 
ter plots of |A0/AF| vs V for different values of 02 . The 
uncertainty is the standard deviation calculated using a boot¬ 
strap method. 


been seen in other numerical experiment^^^ as well as 
in experiments^. To analyse the late stage coarsening 
dynamics, we consider a simple model where the grain 
growth is mediated purely by grain coalescence. The 
grain coalescence follows from the continuous rotation 
of predominantly small grains, which tend to align their 
lattice with larger neighboring grains. Following^U we in¬ 
troduce a characteristic time ti it takes for two grains to 
coalesce. The dynamics of the number of grains can then 
be described by 


I dN _ 1 

N dt ti 


( 4 ) 


Since coalescence between two grains happen when their 
lattices are rotated to be aligned, we assume U to be 
proportional to the inverse of the mean grain rotation. 
Using the conservation equation between the total 
number of grains and mean volume N{t){V) = Vgys and 
differentiating with respect to time and using the coales¬ 
cence assumption, we establish the following relation 


_ 1 dN{t) _ 

{V) dt N{t) dt ^ ' 


( 5 ) 


Assuming (Ad) = Ct~< we find the two solutions 

(m = K)(^) 7 = -l (6) 

(U(t)) = Uoexp 7^-1 

( 7 ) 

From these solutions, we see that if the degree of rota¬ 
tion in the system does not fall off quickly enough grain 
coalescence will lead to exponential coarsening in time. 
In Fig. we observe that for the late stages of the grain 
growth, the decay of (Ad) is in general much slower than 
t~^ and consequently, we expect from Eq. § an expo¬ 
nential grain growth. The fact that we indeed observe 
an exponential grain growth in the simulations supports 
the hypothesis that the grain rotation is the primary me¬ 
diator of grain growth. From Fig. we find that in the 
coarsening regime where 7^0 the proportionality con¬ 
stant C increases for increasing 02 in agreement with the 
inset in Figj^ 

We can further estimate the coupling between grain 
growth and rotation. 


Ad 

AU 


v-P 


( 8 ) 


by assuming that the grain rotation is predominantly 
given by the normal motion of the grain boundary 

/ N dO 

( 9 ) 

The conservation of dislocations then implies thafP^ 
r{t)9{t) = const, or equivalently, 9{t) ~ from 

which it follows that /3 = 4/3 in reasonably good agree¬ 
ment with our numerical results. The fact that simula¬ 
tions agree with this simple estimate might be an indi¬ 
cation that dislocations are indeed conserved up to the 
point where grains coalesce. 

From simulations of grain growth in polycrystalline 
materials, we have identified two distinct dynamical 
regimes. One regime for deep quenching parameters 
where grain rotation quickly is suppressed and therefore 
leads to an overall stagnation of growth. Another regime 
is observed for more shallow quenches where grains con¬ 
tinue to rotate and therefore will be able to align their 
lattices. The alignment eventually leads to the coales¬ 
cence of neighbouring grains and simultaneously allows 
for a few abnormally large grains to form. The grain ro¬ 
tation and coalescence cause an exponential increase in 
grain growth with time. 
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